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Abstract. It is well-known that the free electron Kondo problem can be described 
by a one-dimensional (ID) model because only the s-wave part of the electronic wave- 
function is affected by the Kondo coupling. Moreover, since only the spin degrees of 
freedom are involved in the Kondo interaction, and due to spin-charge separation in 
ID, the universal low energy long distance physics of the Kondo model also arises when 
a magnetic impurity is coupled to the end of a gap-less antiferromagnetic Ji — J2 spin-i 
chain, where Ji{J2) is the (next-)nearcst neighbor coupling. Experimental realizations 
of such spin chain models are possible and using various analytical and numerical 
techniques, we present a detailed and quantitative comparison between the usual free 
electron Kondo model and such spin chain versions of the Kondo problem. For the 
gap-less Ji — J2 spin chain two cases are studied, with zero next nearest neighbor 
coupling, J2 = 0, and with a critical second neighbor coupling, J2 ~ We first focus 
on the spin chain impurity model with a critical second neighbor antiferromagnetic 
exchange J| ~ 0.2412 where a bulk marginal coupling present in the spin chain model 
for J2 < J| vanishes. There, the usual Kondo physics is recovered in the spin chain 
model in the low energy regime (up to negligible corrections, dropping as powers of 
inverse length or energy). At J| the spin chain model is not exactly solvable and 
we demonstrate the equivalence to the Kondo problem by comparing Density Matrix 
Renormalization Group calculations on the frustrated spin chain model with exact 
Bethe Ansatz calculations of the electronic Kondo problem. We then analyze the 
nearest-neighbor model (J2 — 0) where a new kind of Kondo effect occurs due to 
the presence of the bulk marginal coupling. This marginal coupling alters slightly 
the /3-function for the Kondo coupling leading to a slower variation of the Kondo 
temperature Tk with the bare Kondo coupling. In the exact Bethe ansatz solution 
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of this spin chain impurity model ( J2 = 0) Frahm and Zvyagin noted this relation as 
well as the connection to the Kondo problem. Here, by numerically solving the Bethe 
ansatz equations we provide further evidence for the connection to Kondo physics and 
in addition we present low temperature Quantum Monte Carlo results for the impurity 
susceptibility that further support this connection. 



PACS numbers: 75.10.Pq,75.20.Hr,75.40.Mg,03.70.+k,04.20.Jb 
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1. Introduction 

The usual Kondo Hamiltonian [1, 2] contains a Heisenberg interaction between an 
impurity spin and otherwise non-interacting electrons. A simple model takes a free 
electron dispersion relation and a 5 -function Kondo interaction: 

H= f d\[ij\-V^l2m)ij + JK6\^)ij\a/2),lj ■ S^.^J. (1.1) 



Upon expanding in harmonics, this can generally be reduced to a one-dimensional model 
of electrons on a semi-infinite line interacting with the impurity spin at the origin [3]. 
We shall refer to this model as the standard free electron Kondo model (FEKM). A 
closely related model consists of an open gap- less spin-^ Heisenberg chain, defined on 
a semi-infinite line, with one coupling at the end of the chain weaker than the others, 
described by the following Hamiltonian (see Fig. 1): 

L-l 



i=l 



Figure 1. Schematic picture for an open Heisenberg chain coupled to a spin impurity 
(arrow) at the left boundary via a weak antifcrromagnetic exchange J'j^ (dashed bond). 

More generally, it is of interest to add a frustrating second neighbor interaction, J2, 
to this gap-less spin Hamiltonian while keeping the open boundary conditions. In this 
case we arrive at a slightly generalized version of the SCKM: 

L-l L-2 

Ti = Ji Si ■ Si^i + J2 Si ■ Si+2 + 'Himp, 

1=1 1=2 



T^imp — Jx^imp ■ yJiSi + J2S2J ■ (1-3) 
In the following we shall always set Ji = 1. We refer to both spin chain models. 



^ J2J'k 

Figure 2. Schematic picture for the open frustrated Heisenberg chain coupled to a 
spin impurity (arrow) [model (1.3)]. Thin lines show the fr'ustrating second neighbor 
couplings J2. 

Eqs. (1.2) and (1.3), as the Spin Chain Kondo Model (SCKM). As J2 is increased a 
critical point J| ~ 0.2412 [4] is reached beyond which the spin chain enters a dimerized 
phase [5] with a gap and the relation between Eq. (1.3) and Kondo physics no longer 
holds. 
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The purpose of this paper is to explore in some detail the relationship between 
the FEKM and the SCKM. Such a relationship was pointed out in [6, 7] using field 
theoretical techniques and numerical solutions of finite chains. In [8] , hereafter referred 
to as FZ, it was pointed out that the spin Hamiltonian of Eq. (1.2) (or Eq. (1.3) with 
J2 = 0) is Bethe ansatz integrable and a number of its properties were determined 
exactly, establishing some connection with the usual (free electron) Kondo model. 
However, certain clear differences between the two models are evident from the exact 
solution. In particular, the characteristic energy scale, the Kondo temperature, behaves 
exponentially at small Kondo coupling, T/^ oc e~'^l'^*^ (where c is a constant) in the 
FEKM but has an unusual exponential square root behavior [8] in the SCKM at J2 = 0, 
Tft; oc e~^l^^ . This behavior is likely generic for all J2 < J|. However, as we shall 
see, once we tune J2 to the special critical point J2 ~ 0.2412 in the SCKM the usual 
exponential behavior Tk oc e~'^l'^'K is recovered. 

Here we will explain the subtle differences in the low energy properties of the FEKM 
and SCKM using a renormalization group/field theory approach. The differences arise 
from a marginally irrelevant interaction in the spin chain model, with marginal coupling 
constant g. This is a hulk interaction, which does not occur in the FEKM and is unrelated 
to the Kondo interaction. At ^^e marginal coupling constant vanishes, (7 = 0, and 
the usual Kondo physics is exactly recovered in the SCKM. 

In the strong coupling limit of the usual Kondo model, the impurity spin is 
'screened' in the sense that it will form a non-magnetic singlet with the rest of the 
system. The screening is expected to occur on a characteristic length scale, defining the 
so called screening cloud: 

ei, = t;/r^cxei/"-. (1.4) 

Here v is the characteristic velocity of the low energy excitations. Physical quantities 
with zero anomalous dimension are expected to show scaling with I L. Strong evidence 
for the occurrence of Kondo physics in the low-energy sector of the SCKM can therefore 
be established if estimates, not only of (or equivalently Tk), but of the entire scaling 
functions coincide in the SCKM and FEKM. We show that this is clearly the case for 
the SCKM at J2 = J| while for the SCKM at J2 = the logarithmic corrections arising 
from the marginal coupling constant g leads to scaling violations and as mentioned 
above, the divergence of Tk and therefore ^k acquires an unusual exponential square 
root behavior. Still, the strong coupling limit remains that of a fully screened impurity 
for J2 = (^ > 0). 

Our conclusions will be confirmed in detail through large-scale numerical 
calculations. To this end, we have made extensive use of the two known exact Bethe 
ansatz solutions, for the FEKM [9, 10] and for the SCKM at J2 = [8]. Numerical 
results are then obtained by solving the two sets of Bethe ansatz equations (BAE) for 
finite systems of size L at T = 0. At finite T, we also obtain quantum Monte Carlo 
(QMC) results for the SCKM at J2 = 0. For general J2, and in particular at J| where 
the relation between the SCKM and FEKM is most direct, the SCKM is not exactly 
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solvable and a sign problem appears in QMC simulations due to the frustrating term 
J2. Fortunately, high precision results for the SCKM at T = in this regime (J2 7^ 0) 
can be obtained using Density Matrix Renormalization Group (DMRG) methods. 

These results are useful in two ways. First of all, the SCKM has distinct advantages 
over the free electron version for numerical simulation since the Hilbert space grows much 
more slowly with system size, L, as 2^ rather than the 4^ behavior of the free electron 
model. (On the other hand, the Numerical Renormalization Group technique provides a 
powerful alternative numerical approach for the free electron model.) In particular, we 
presented extensive DMRG results on entanglement entropy for the Kondo model using 
the spin chain version in [11, 12, 3]. In addition, experimental realizations of the SCKM 
are possible, in particular in muon spin resonance experiments [13]. In experiments on 
spin chains the coupling J2 usually cannot be tuned and a detailed understanding of the 
relation to the FEKM as J2 is varied over experimentally realistic values is therefore of 
considerable interest. 

This paper is organized as follows. In section 2 we briefly review the field theory 
and renormalization group (RG) treatment of both free electron and spin chain Kondo 
models. We then study the interplay of the marginal bulk interaction and Kondo 
interaction in the SCKM, explaining the origin of the exponential square root form 
of the Kondo temperature using a renormalization group approach. In section 3 we 
present results on the SCKM with a (bulk) next nearest neighbor interaction, Eq. (1.3), 
tuned to the critical point, J| ~ 0.2412, separating gap-less and dimerized phases. The 
presence of the frustrating coupling, J2, allows the bulk marginal coupling, g to be 
tuned and at the critical point, J|, it vanishes, = 0, where excellent agreement is 
obtained with the FEKM. A precise relationship between Jj^ in the SCKM and the 
Kondo coupling, A/^, in the FEKM is obtained in this case which is needed for detailed 
comparisons [3, 12]. This relation is established in sub-section 3.1. In sub-section 3.2 
we then turn to a discussion of the finite-size corrections to the triplet-singlet gap in 
the FEKM and SCKM (at Jif) and the resulting scaling with ^k/L. This allows for 
a convenient way of extracting and we show that estimates of and the entire 
scaling function coincide in the FEKM and SCKM at J2 = J|. We end section 3 by 
deriving the finite-size corrections to the ground-state energy in sub-section 3.3. As 
previously shown [3, 12], a universal term proportional to the Kondo length scale can 
be identified in the ground-state energy at intermediate values of the Kondo coupling. 
This allows for an independent and absolute estimate of in the SCKM in excellent 
agreement with those obtained from the scaling of the singlet-triplet gap. In section 4 
we turn to a discussion of our results for the nearest neighbor SCKM, Eq. (1.2), with 
J2 = where (7 > and logarithmic corrections arising from g are present. In sub- 
section 4.1 we first calculate the impurity susceptibility of the SCKM (J2 = 0) at finite 
temperature using QMC and compare it to the known result for the FEKM obtained 
from the BAE. In sub-section 4.2 we discuss the singlet-triplet gap for the finite size 
system. The result is compared with weak and strong Kondo coupling perturbative 
results, supplemented by weak couphng perturbation theory in g for the SCKM. This 
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is shown to be a convenient quantity for extracting the Kondo screening cloud size, 
^K- It compares well to the expected fs/T^- where is the characteristic velocity. 
We then discuss the finite size corrections to the ground-state energy for the SCKM 
with finite and (7 > in sub-section 4.4. We close section 4 by discussing Bethe 
ansatz results of the T = magnetization [8], M{H), in sub-section 4.5. We observe 
that, in addition to the standard scaling function of H/Tk which occurs in the FEKM, 
additional scaling violating terms appear which can be expanded in g{H). Finally, in 
section 5 we summarize our results by showing Tk{J'k) as obtained for the SCKM both 
at J2 = and J2, and compare to the formula obtained by FZ and to that obtained 
from RG calculations in section 2. 



2. The Kondo effect: A one dimensional problem 

2.1. Low energy theory 

Due to the spherically symmetric free particle dispersion relation of the 3D model in 
Eq. (1.1), we may expand the electron annihilation operators in spherical harmonics. 
Due to the 5-function form of the Kondo interaction, only the s-wave harmonic interacts 
with the impurity. (See, for example, [14] for more details.) Assuming a weak Kondo 
coupling, we may linearize the dispersion relation near the Fermi momentum, kp, 
yielding the ID low energy effective Hamiltonian: 



H='^ I dx 



27r Jo 



00 



(2.1) 

Here the left and right moving Fermi fields are functions of {vpt + x) and {vpt — x) 
respectively and obey the (unconventionally normalized) equal time anti-commutation 
relations: 

{^L/R{x),^lj,iy)} = 2tx5{x - y), (2.2) 

and the boundary condition at the origin: 

^|JL{x = 0) = -^Ijr{x = 0). (2.3) 

[We have adopted the normalization conventions for the fermion fields and the definition 
of \k used in [14]. \k = vJk where u is the density of states (per spin) u = kprn/ (27r^).] 
A convenient way of analyzing the ID model is to bosonize, since this separates 
spin and charge degrees of freedom. First we take advantage of the boundary conditions 
(2.3) and the fact that 'ipL/Rit^x) = ipL/Ri'Vpt ± x), to write the right movers as the 
analytic continuation of the left movers to the negative axis: 

M-x) = -^Rix), > 0). (2.4) 
Then the Hamiltonian, written in terms of left-movers only, becomes: 

H='^ r dxti^l^ijL + vfXk4^1{Q)^M^) ■ S. (2.5) 
2vr dx 2 
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It is possible to write the kinetic energy in a form quadratic in charge and spin currents: 



i,«t/^N„/, ( ^\ 7 „/,"t a 



(Here repeated indices are summed.) The Hamihonian becomes: 

vpXKJm-S. (2.7) 



H = — I dx 

2ti 



oo 



oo 



1 1 - - 

- : JlJl{x) : ■ Jl 



The spin and charge currents operators may be expressed in terms of spin and charge 
boson fields and only the spin bosons couple to the impurity, S. The kinetic energy 
terms just correspond to free massless (charge and spin) boson kinetic energies. 

2.2. Renormalization group 

For a weak bare coupling, we can estimate S,k as the length scale at which the effective 
coupling becomes 0(1). From the weak coupling /3-function, the variation of the effective 
Kondo coupling as we vary the length scale is: 



Here Xk{L) is the effective coupling at length scale L. Alternatively, we may identify 
L with Vp/E where vp is the velocity and E is a characteristic energy scale such as a 
temperature or magnetic field. 

We may integrate this equation from a length scale of the order a lattice constant, 
a, up to with the effective Kondo coupling varying from its bare value A^^- to a value 
c, of 0(1). Assuming A^^ <C 1, this gives approximately: 

^ [ae-'/^y^ -^exp[l/A?,]. (2.9) 



Note that the first factor, in square brackets is simply a constant. Including higher 
order terms in the /3-function only leads to relatively small corrections of the form: 

= exp[l/A;,] [1 + 0{Xk)] . (2.10) 

(Here we have dropped the subscript from Xk which refers to the bare value of the 
Kondo coupling.) We may solve for the effective coupling at scale L in terms of at 
weak coupling (i.e. L -C ^k)'- 

2.3. Fermi liquid theory 

At low energies and long length scales, the effective Kondo coupling becomes large, and 
the effective Hamiltonian flows to the strong coupling fixed point. We may think of 
this fixed point as one where the impurity spin is "screened", i.e. it forms a singlet 
with a conduction electron. The remaining electrons behave, at low energies and long 
length scales, as if they were non-interacting, except that they obey a modified boundary 
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condition reflecting the fact tliat tliey cannot break up tlie singlet by occupying the same 
orbital as the screening electron. This modified boundary condition corresponds to a 
tt/2 phase shift. Correspondingly in the spin chain Kondo model, the impurity spin gets 
"adsorbed into the chain" and no longer behaves like a paramagnetic spin at low energies 
and long distances. The leading corrections to this low energy long distance picture are 
described by lowest order perturbation theory in the leading irrelevant operator at the 
strong coupling fixed point. This is an interaction between the remaining conduction 
electrons, near the screened impurity. (It doesn't involve the impurity itself since it is 
screened and doesn't appear in the low energy effective Hamiltonian.) 

This leading irrelevant operator is Jl ■ Jl(0)[15, 16]. From (2.7), we see that this 
is proportional to the spin part of the free electron energy density, 7is,L(0). The energy 
density has dimensions of (energy) /(length) so the corresponding coupling constant in 
the effective Hamiltonian must have dimensions of length. On general scaling grounds 
we expect it to be proportional to ^k- The precise constant of proportionality simply 
corresponds to giving a precise definition of what we mean by ^k- We adopt the 
convention: 

H^nt = -{7l^K)n,,L{0)- (2-12) 

Here the subscripts s and L are a reminder that this is the spin only part of the energy 
density for left movers. Note that if we start with a system of length L imposed OBC 
(with left and right movers), then we can map into a system of periodic length 2L 
with left movers only. For the purpose of doing first order perturbation theory in Hint 
for quantities like the susceptibility, specific heat or ground state energy, which are 
translationally invariant in 0*^ order, we may replace [14] Hint by: 

H,nt^-^ j'^HsA^)- (2.13) 

This is equivalent to a length dependent reduction of the velocity: 

VF^VF 1-^ . (2.14) 

This then implies that the T = susceptibility, which is L/ {2ttvf) in the absence of the 
Kondo impurity becomes: 

L L ex 



X 



27rt;^[l - 7r{i^/(2L)] 2nvF ^Vp 

L 1 



2nv ■ 4T, P^^^' 

Thus the zero temperature impurity susceptibility is 1/{ATk)- It is this form of the 
impurity susceptibility, simply related to the high temperature, free spin behavior, 
1/(4T), which motivates the definition of (and hence Tk = vf/^k) implied by 
(2.12). We note that this interaction Hint is present even in the absence of an impurity, 
for free fermions but then the coupling constant is of order a lattice constant. Similarly, 
it is also present for the spin chain with no impurity (i.e. J'j^ = 1) with a coupling 
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constant of order a lattice constant. The effect of a weak Kondo coupling is to make 
this coupling constant large. We emphasize that this precise choice of definition of 
has no physical consequences. The power of Fermi liquid theory is to predict not only 
the form of low energy quantities but also ratios of various low energy quantities such as 
impurity susceptibility, impurity specific heat, resistivity, etc., corresponding to various 
generalized Wilson ratios. 

2.4- Field theory approach to the Spin Chain Kondo Model 

A field theory approach to the spin chain Hamiltonian of Eq. (1-2) is obtained by 
bosonization. The spin operators have uniform and staggered parts: 

Sj [Jiiaj) + jR{aj)]/{2Tx) + (-Inconstant ■ n(aj), (2.16) 

where a is the lattice constant and Jl/Ri n vary slowly on that scale. The left 
and right moving current operators are equivalent to the corresponding operators in 
the free fermion model. On the other hand, the staggered spin density, n, with a 
scaling dimension of 1/2, has no counterpart in that model. The low energy effective 
Hamiltonian is simply that of the spin part of the free electron Hamiltonian, 

Ho = ^ + Jl]. (2.17) 

67r Jo 

Here Vs is the spin velocity. In addition, there are various irrelevant operators appearing 
in the effective Hamiltonian. The only one which is important at long distances and low 
energies is a marginally irrelevant interaction, 

H,r.t = -^-^ r Jl-Jr. (2.18) 
27r Jo 

The bare coupling constant, go, has a positive value of 0(1) but renormalizes to zero 
at low energies and long lengths scales. As a frustrating next neighbor interaction, J2 
is turned on go decreases, passing through zero at a critical coupling, J2 ~ 0.2412. For 
still larger J2, go becomes negative, and hence marginally relevant. 
The RG equation for the bulk marginal coupling constant, g, is: 

+ (2.9) 

Integrating from a scale, Lo of order a lattice constant where the bare coupling has a 
value > to an arbitrary scale L gives: 

'^'■^^ da 

-^=ln(L/Lo). (2.20) 

It is known that g{L) flows, at large length scales, to the weak coupling region where 
the expansion of Eq. (2.19) becomes valid. For large L, the integral in Eq. (2.20) is 
dominated by the small g region giving: 
1 1 

- In g{L) = ln(L/Lo) + constant = ln(L/Li). (2.21) 



90 



9{L) 2 
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The constant Li cannot, in general, be determined by elementary means. However, if 
the bare coupling, Qq is weak then: 

ln(L/Li) ^ ln(L/Lo) + -, (2.22) 

^0 

so Li ^ Lq. For the nearest neighbor Heisenberg model, Li is 0(1) and can be 
determined by comparing to the Bethe ansatz solution (for Jj^ = 1) [17]. In this case, 
the best choice of Li depends on precisely what quantity is being calculated. In any 
event, at sufficiently large L, where L ^ Li, we may approximately solve Eq. (2.21) to 
obtain: 

,(L).-i---l^-N^/^-)l (2.23) 
' ln(L/Li) 2 \n\L/Li) ^ ' 

Note that the /3 functions for the Kondo coupling, \k{L) and bulk marginal coupling, 
g{L) are the same, to cubic order except for the sign of the quadratic term. Hence the 
expressions, Eq. (2.11) and (2.23) look similar. There are crucial differences however. 
Eq. (2.11) is only valid at relatively short distances (high energies) where L -C ^k'-, 
\k{L) grows with increasing L. Eq. (2.23) is only valid at long distances (low energies), 
L ^ Lq for the nearest neighbor Heisenberg model; g{L) decreases with increasing L. 

For go < 0, there is a gap of order \/\go\e~^^^^°^, corresponding to the spontaneously 
dimerized phase. The critical point, where go = is thus an especially good point for 
numerical studies since, for g^ > the slowly decreasing marginal coupling constant 
complicates finite size scaling. 

Weakly coupling an impurity spin to a single spin at a generic point in a spin chain 
yields a model [6, 7] with essentially no relation to the standard Kondo model. On the 
other hand, the model of Eq. (1.2) with the impurity spin coupled at the end of the 
chain, is quite different. To see this, first consider the limit J'j^ = g = 0. The open 
boundary condition leads to: 

4(0) = /L(0)cxn(0). (2.24) 

Thus the SCKM is equivalent to the FEKM for some non-trivial constant of 
proportionality between Kondo couplings, and with the Fermi velocity replaced by the 
spin velocity, ignoring the irrelevant operator. 

The low energy effective Hamitonian of the SCKM Eq. (1.2) has three main terms: 



if = — / dx- 
2n ./_L 3 



L 



Jdx) 



gvs 

27r 



L 



dxJiix) ■ Jr{x) + XRVsJiifi) ■ S. 

(2.25) 

The first term, which we denote Hq, corresponds to a free boson Hamiltonian, the 
second term is the marginal bulk interaction, and the third term stands for the Kondo 
like interaction at the boundary with \k oc J'j^. Using the boundary condition, 
JL{t,x = 0) = JR{t,x = 0) we conclude that we can regard Jr{x) as the analytic 
continuation of Jl{x) to the negative x-axis: 

Jnix) = M-x). (2.26) 
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Hence we may write the Hamiltonian as: 



H = Ho-^l dxJL{x)-JL{-x) + XKVsJL{0)-S. (2.27) 



2vr JO 

When g = this is the same field theory that describes the FEKM, Eq. (2.7), with vp 
replaced by Vg, the spin velocity. A non-zero g leads to some differences, but since g is 
marginally irrelevant, its effects becomes progressively less important at long distances 
and low energies. 

We now consider the renormalization group equations for the model with both g and 
Xk 7^ 0. The presence of the Kondo interaction, which is only at the boundary, cannot 
change the renormalization of the bulk interaction, g. On the other hand, the reverse is 
not true; the bulk marginal interaction has an important effect on the renormalization 
of the Kondo coupling constant. This can be viewed as resulting from the fact that 
the boundary operator Jl(0), appearing in the Kondo interaction, which has a scaling 
dimension of 1 when g = picks up an anomalous dimension of first order in g. This 
anomalous dimension was calculated in [18]: 

x=l-g + ... (2.28) 

[Note that, in [18], g was normalized differently, so that gAQ = (v^/47r)(7.] This 
observation determines the RG equations: 

i^^«A, + Ai..... (2.29) 

To solve (2.29) we may substitute g{L) from (2.23), yielding: 
dXK Xk ,2 



This differential equation can be solved exactly by defining a new effective coupling 
constant: 

^'^) - m ''''' 

which obeys the RG equation: 

^ = Pl„(L/LO. (2.32) 

Setting Xk{^k) = 1, and using Eq. (2.23), one can obtain the renormalized expression 
for the effective Kondo coupling 

' (1/2) \n\U/Li) - (1/2) \n\L/L,) + HU/U) ' ^ ' ' 

This gives, for the Kondo length scale ^k, the unusual dependence on the bare Kondo 
coupling X^j^: 



= Li exp 



c+J2\n{Lo/L,)/X% + \n\Lo/U)+c^ 



(2.34) 



where c is a positive constant of 0(1). 
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Clearly, Xk{L) is not a scaling function of ^k/L. Nevertheless, when the bare 
marginal coupling qq is very small, which is true when J2 is tuned close to the critical 
point, the length scale 

Li ~ exp(-l/(?o) < Lo (2.35) 

and then 

Assuming A^^- -C l/ln(Lo/Li), with a fixed bare marginal coupling Eq. (2.34) 
may be approximated: 



U ^ Li exp ^2 ln{Lo/L,)/X%. (2.37) 

Note that this dependence of on the Kondo coupling, oc exp[b/ ^/)^] (where b 
is a constant), gives a much shorter Kondo length scale than in the usual Kondo effect 
which occurs when Qq = 0: oc exp[l/A5^-]. Also note, from Eq. (2.35) that the limit 
of zero marginal coupling corresponds to Li 0. In this limit, ln(a/Li) ^ l/A^^- so we 
may Taylor expand the square root in (2.34) giving ~ expfl/A^s;]. However, (2.34) 
and (2.37) only hold when ^ Li where Li is roughly the scale at which the effective 
marginal bulk coupling, g{L) starts to become small. 

We expect physical quantities with zero anomalous dimension to be generalized 
scaling functions. This means that they can be written as functions of ^k/L and g{L) 
only. In general, they should have a Taylor expansion in the renormalized coupling g{L): 

00 

f{XK,9,L) = Y,9^{L)fn{^K/L). (2.38) 

n=0 

The leading term, fo{C,K/L) should be the same as in the model with g = 0, the FEKM. 
Note that, due to the presence of the g"'{L) ~ 1/ ln"(L), / is not a pure scaling function 
in the sense that it does not depend only on ^k/L, but since g{L) ^ as L ^ 00, 
^k/L scaling becomes a better and better approximation at large L. More generally, 
the inverse length 1/L in Eq. (2.38) may be replaced by an energy scale such as a 
magnetic field. 

The modified (exponential square root) expression for the dependence of the Kondo 
temperature on Kondo coupling and the slowly decreasing correction, gfi, (7^/2, . . . are 
the two main differences between the FEKM and the SCKM. In the following sections we 
demonstrate the correctness of this assertion by considering various physical quantities. 



3. Kondo effect in the SCKM at 

We first discuss the frustrated spin chain model (depicted in Fig. 2) with the Hamiltonian 
Eq (1.3). As we have already stressed, in this case the marginal coupling is exactly 
zero, g = 0, and the mapping to the FEKM becomes exact, up to strictly irrelevant 
interactions. Hence, in this case, by determining various constants numerically, it is 
possible to establish a precise mapping between the low-energy, long-distance sectors of 
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the FEKM and SCKM. In sub-section 3.1 we numerically determine Vg in the SCKM 
as well as the constant of proportionality between in the SCKM and \k = vJk in 
the FEKM. Then, in sub-section 3.2, we directly compare Bethe ansatz results for the 
FEKM with DMRG results for the SCKM at J|. Finally, in sub-section 3.3 we show how 
at intermediate couplings, Jj^, for the SCKM can be determined from the finite-size 
scahng of the ground-state energy in the SCKM. 

First consider the model with J'j^ = 0. At a general site, j, the low energy degrees 
of freedom of the spin operators are represented as: 

~ ^[JLia-j) + JRia-j)] + (-l)-'constant x n{aj). (3.1) 

Here a is the lattice spacing and n{aj), the alternating part of the spin operators can be 
written in terms of the spin boson field in a non-linear way. At the end of the chain, due 
to the open boundary condition, we find that n{x) —>■ constant x Jl(0) and therefore: 

S2 + J2S3^CJl{0), (3.2) 

where C is a non-universal constant, depending on the second neighbor coupling, J2 in 
the Hamiltonian. [C has dimensions of inverse length, and is proportional to the inverse 
lattice spacing.) Now including a weak coupling, J'^ to the first spin, the low energy 
effective Hamiltonian becomes the usual FEKM with the replacements: 

vp Vs, vpXx CJ'j^. (3.3) 
3.1. Numerical parameters for the spin chain problem 

Velocity of excitations: In order to make the above correspondence quantitative, one 
needs to determine exactly the spin velocity Vs at the critical point J2 = J|. We can use 
the fact that the finite size scaling of the ground state energy for the uniform SCKM 
with = 1 is simply given by [19] 

EoiL) = eoL + - ^ + g + g + 0(1/1% (3-4) 

eo, ei and 02, 03 are non-universal numbers, whereas the 1/L term is universal, 
proportional to the central charge (c = 1 here) and the spin velocity Vs- Note that, 
Vs will depend on J2 and even though it is known analytically from the Bethe ansatz 
solution available at J2 = we need to determine it at J| where no exact results are 
available. DMRG results for Eo{L), computed on chains up to L = 1000 sites, are 
shown in Fig. 3 and fitted to Eq. (3.4). We have checked that fitting with a M^'^-order 
polynomial form in 1/L does not change the estimate of Vg whenever M > 3. Thus we 
found for the spin velocity 

Vs = 1.174(1). (3.5) 

This is in good agreement with previous estimates [6]. 

Kondo coupling: Of course, to predict for the spin chain model we need to 
exploit Eq. (3.3) to relate to Jj^, the weak coupling of the spin at the end of the 
chain, where the constant C is determined by (3.2). We emphasize that this constant C 
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Figure 3. Finite size scaling of the ground state energy Eq{L) of the SCKM at the 
critical second neighbor coupling J| = 0.2412 shown for chains up to L = 1000 sites. 
DMRG results obtained keeping m ~ 512 states (open symbols) are compared to the 
finite size scaling form (3.4) up to third order (full line). In the inset, Vs obtained 
from polynomial fits of the form Eq{L) = EqL + ei — ttVs 

/{24:L) + J2iL2 ai/L' is plotted 

versus M. 



can be determined in the theory with zero Kondo couphng, = and we can therefore 
study correlations in chains of even length. A convenient way of determining C is to 
measure a long-distance correlation function in the spin chain and compare it to the 
corresponding correlation function in the continuum field theory. The simplest choice 
appears to be the end-to-end equal time correlation function: 

((5^2 + J2S,) ■ {Sl + J2Sl^i)) ~ C^iMO) ■ ML)). (3.6) 
The current correlation function is simply that of the free fermion model with the current 
defined in Eq. (2.6) and the fermion fields normalized as in Eq. (2.2). For an infinite 
chain this implies: 

{Ur) ■ Mr')) = ^ (3.7) 

However, we need the correlation function for a finite strip. Using the fact that the open 
boundary conditions on a finite strip of length L are equivalent to periodic boundary 
conditions for left movers only on a circle of circumference 2L, we obtain: 
^ ^ 3 

(8L^/7r^) sm [7r(r — r) / 1L\ 
giving the end-to-end correlation function: 

37r2 

{J,{Q) ■ J,L)) = - — . (3.9) 

As shown in Fig. 4 this correlation function can be determined very accurately using 
DMRG methods, from which we determine: 

{02 + J2S,) ■ {Sl + J2SlM) - -^^2^11^. (3.10) 




Figure 4. DMRG results keeping m = 256 states for {{SI + J2SI) ■ {SI + J2Sl_{)) for 
the SCKM at J| (open circles). Results are shown for L = 18, 20, 30 . . . 110. The sohd 
circles indicates a polynomial fit with the leading term of the form — 0.891462L~^. 



From Eqs. (3.6), (3.9) and (3.10) we obtain: 

0.891462 = ^— , C = 0.850054 (3.11) 

8 

From (3.3), using Vg ~ 1.174(1), we finally obtain the proportionality constant between 
Kondo couplings in spin chain and fermion models, Jj^ and A^: 

J'k = {vs/C)\k = 1.38 Xk. (3.12) 

3.2. The triplet- singlet gap and Bethe ansatz results for the FEKM versus DMRG 
results on the SCKM at J2 

A convenient diagnostic for testing the correspondence between the FEKM and the 
SCKM at J2 is the finite size behavior of the gap between the singlet ground state and 
lowest triplet exited state (for even L). Let us first consider the models with no Kondo 
coupling. Thus we consider either a free fermion chain or a Heisenberg spin chain with 
open boundary conditions. Noting that the allowed wave- vectors are k ^ nn/L, we see 
that the excitation energy of a triplet particle-hole excitation in the free fermion chain 
is: 

AsT = (3.13) 

The same result holds for the spin chain with open boundary conditions [6]. The 
extension of this result for ^ L, is easily obtained from Fermi liquid theory, by 
the usual replacement Eq. (2.14): 
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Figure 5. Scaling function Q{L/^k) (see Eq. (3.16)) obtained from tlie Bethe Ansatz 
solution of the one-channel FEKM Hamiltonian [1] with various values for the Kondo 
interaction strength Jk, as indicated on the plot. The scaling plot is constructed using 
the strong coupling limit curve [long-dashcd green line Eq. (3.14)] as a support for the 
rest of the collapse. The spin chain model (1.3) at J| results (star) arc adjusted to 
collapse onto the universal curve. Left inset: in the weak coupling regime, the ratio 
O ~ l/ln(5if/L), as show by the linear fitting dotted line. Right inset: in the strong 
coupling regime, the FLT result 6 = 1 — describes perfectly the behavior (long 
dashed hue). 



In the weak coupling limit we may calculate Ast in perturbation theory in the Kondo 
coupling, \k- We consider a spin chain with an even number of sites, or equivalently 
a free fermion chain with an odd number of fermions, coupled to the impurity spin so 
that the ground state is a spin singlet. In the limit of zero Kondo coupling the ground 
state of the non-interacting chain is a spin doublet. First order perturbation theory in 
Xk couples this doublet to the impurity spin with a coupling constant nVs/L, giving a 
singlet-triplet splitting of /S.st ~ {j^Vg/ L)\k- We expect that, as usual, higher order 
corrections will replace A^^ by its renormalized value at scale L, Eq. (2.11), giving: 

TTf, 1 

- XMWI) ■ ("^' 

Hence, in general we expect the quantity 

e = — A5T, (3.16) 

TlVs 

to be a universal scaling function of L/S,k- The same function should occur for the 
FEKM and for the SCKM up to higher order corrections from irrelevant operators. 
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J'k 


0.6 0.45 0.3 0.25 0.2 




5 17 130 400 1800 



Table 1. Kondo length scale estimates for the SCKM at J|, Eq. (1.3). 

Introduced 25 years ago, independently by Andrei [9] and Wiegman [10], the exact 
diagonalization of the Kondo Hamiltonian (2.5) based on the Bethe Ansatz provides a 
powerful method to access several physical observables like the uniform susceptibility, 
the field induced magnetization, the specific heat, the thermodynamic entropy. Similarly 
to what was previously done using the Bethe Ansatz solution of the unfrustrated chain, 
one can add an external magnetic field along the 2;-axis B = Be^, which couples to 
the total spin operator S^^^. Then one can solve numerically for finite length chains 
the Bethe Ansatz equations [9, 10] in any S^^^ sector and compute the singlet triplet 
excitation gap. 

We have calculated Q{L/^k) for the FEKM by solving the set of coupled Bethe 
Ansatz equations for various values of the Kondo coupling J, where the convention of 
Ref. [1] was followed: J = |A/<. In Fig. 5 are shown the results for the ratio Q{L/^k), 
successfully compared with the strong and weak coupling predictions Eqs. (3.14) and 
(3.15). 

We have also compared these results to DMRG calculations performed on the Kondo 
spin chain model with J2 = J|, (1.3) for J'j^ = 0.2, 0.25, 0.3, 0.45, 0.6. The ratio 9 
is determined for L = 100 sites and the DMRG data points (shown as ^) have been 
superimposed on the Bethe Ansatz results on Fig. 5, adapting L/^k{J'k) i'^ order to get 
a good collapse. The resulting values for the Kondo length scales are listed in table 1. 

The fact that the numerical estimates of not only the length scale ^k-, but of the 
entire scaling function coincide in the FEKM and the SCKM at J2 = J2 is very strong 
independent evidence for the occurrence of Kondo physics in the low-energy sector of 
the SCKM at J2 = J|. 

3. 3. Finite size scaling of the ground state energy in the SCKM at J| 

One way of estimating ^k, in the case L, is from a universal correction to the 

ground state energy. This may be calculated from the Fermi liquid theory Hamiltonian 
of (2.12). It is well known that the ground state energy of the open spin chain has the 
form: 

EGs = eo + e,L-^ + ... (3.17) 

where eo and ci are non-universal but the third term is universal, depending only on 
Vs and . . . indicates terms that drop off faster with L. We may include an additional 
boundary term from first order perturbation in the FLT interaction of (2.12) by replacing 
Vg by its shifted value in (2.14), giving: 

6Egs = 7r\sU/{^^L^)- (3.18) 
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We emphasize that such a 1/L^ term is always present, even for the uniform chain with 
open boundary conditions (J^ = 1) but that the factor oi is of order the lattice 
spacing in that case. This term obtains a large coefficient at weak Kondo coupling. 
Combining these contributions we find for the ground-state energy: 

Eos{L^ J'k) = eoiJ'K) + (e^ + ^. (3.19) 

Please note that, while the constants ei,e2 describe bulk behavior and hence are 
independent of Jj^, the surface term cq does depend on J^. Secondly, the difference 
between 62 and ^/^( Jj^ = 1) is largely a matter of convention. Thirdly, this expression for 
the ground-state energy neglects logarithmic corrections arising from the bulk marginal 
operator, g, assumed to be zero since J2 = The more general case with J2 < J2, g > 
will be considered in section 4.4. 

Using DMRG results for the ground-state energy we now attempt to determine 
^k{Jk) using Eq. (3.19). As outlined in Appendix B we can in this case eliminate the 
3 first terms in Eq. (3.19) by using a form of Richardson extrapolation. We then arrive 
at: 

= D4{L) + 0{L~^) (3.20) 

where the coefficient c|(L) is known. Plotting E'^'^Y^i should then yield the desired 
constant, D. Plots of E^^^cKL) versus 1/L are shown in Fig. 6. For = 0.6,0.8, 1.0 

I — I — I — I — I — I — I — I — I — I — 

o j'k=o-3o : 

□ J'^=0.37 

O J'k=0.41 J 

A J'|^=0.45 

<1 J'^=0.525 

T JV=0.60 

p. J\=0.80 

+ J'k="'00 



"0 0.01 0.02 0.03 0.04 0.05 0.06 

1/L 

Figure 6. DMRG results keeping m = 256 states (m = 512 for = 1) for 
E'-^^L, Jk)/4{L) for the SCKM at J| for a range of plotted versus 1/L. For the 
larger the results quickly approach a constant as 1/L 0. Results are obtained 
from the ground-state energy for systems with even size. 

it is easy to extract the constant D using a simple polynomial fit. However, for smaller 
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J'k 


1.000 


0.800 


0.600 


0.525 


0.450 


0.410 


0.370 


0.300 


D 


0.24093 


0.54043 


1.4853 


2.5426 


4.76 


7.14 


11.5 


~ 35 




=0.65 


1.89 


5.81 


10.19 


19.4 


29.3 


47 


~145 



Table 2. The numerically determined values for D for the SCKM at J| and the 
resulting 

we have found it necessary to fit to a (2,2) rational polynomial to improve the 
extrapolation. For Jj^ smaller than J'j^ = 0.3 it is no longer possible to extrapolate using 
rational polynomials and a reliable numerical determination of D becomes impossible. 

As mentioned above, we expect to be C(l) for = 1. In previous work [3, 12] 
we used Fermi liquid theory to determine ^k{J'k = 1) = 0-65 and here we use this result 
to calibrate our data. In particular, we use this result to determine the coefficient 62 
using our previous estimate for the velocity, Vg, finding 62 = 0.0841. The remaining 
are then trivial to obtain. In table 2 we list D{J'j^) along with the resulting ^k- These 
independent (and absolute) estimates of $,k very nicely agree with the ones previously 
obtained in table 1 clearly establishing that the same length scale can be extracted 
from the scaling of the ground-state energy and the singlet-triplet gap as one would 
expect from Kondo physics. 

4. Kondo effect in the nearest neighbor SCKM, J2 = 

We now turn to a discussion of the nearest neighbor SCKM, Eq. (1.2), in the absence of 
any second neighbor coupling (J2 = 0). As outlined, the bulk marginal coupling in the 
SCKM, g, is in this case non-zero, g > 0. An exact Bethe ansatz solution of this model 
was developed by FZ [8]. We start by showing QMC results for this non- frustrated 
model, before presenting BA results for the singlet triplet gap. Then we discuss the 
connections with the usual Kondo physics. 

4.I. Impurity susceptibility 

We now want to compute the spin susceptibility of the SCKM defined in Eq. (1.2). We 
take advantage of the non frustrated nature of this antiferromagnetic spin chain model 
to perform large scale QMC simulations using the Stochastic Series Expansion of the 
partition function in a loop algorithm framework (see Ref. [20] for details). 

4-1.1. Quantum Monte Carlo results While the QMC simulations are done on finite 
size systems, on can still define the impurity susceptibility by 

Xhnp = lim [x{L + 1, J'j,) - xm J'k = 1)] , (4.1) 

where L + 1 corresponds to a system of L spins in the bulk coupled to 1 boundary 
impurity with J'j^, and L corresponds to a pure chain of L spins (see Fig. 1). The QMC 
simulations have been carried out for L = 512 at low temperatures for various Kondo 
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exchange couplings, J^. While we are not in the thermodynamic limit L — * oo, we still 
expect for such a rather big system size that one can capture all essential features of the 
Kondo physics present in the spin chain model. Results for x(513, J'^) are displayed in 
Fig. 7 where a clear upturn is visible at low temperature. 
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Figure 7. Total susceptibility x(T) obtained by Quantum Monte Carlo simulations 
of the SCKM (J2 = 0) with chains of length L ~ 512 coupled to an impurity spin with 
various Kondo couplings indicated as various symbols on the plot. Results for a 
pure open chain are also shown (solid curve). 



As in the usual Kondo problem, one expects the impurity spin to be essentially free 
at very high temperatures T ^ T^-, leading to a Curie-like divergence ~ 1/T. On the 
other hand, when the temperature is decreased, the effective Kondo coupling starts to 
grow and the impurity eventually becomes screened when T ^ Tk-, forming a strongly 
entangled spin singlet with the bulk, and thus resulting in an absence of features in the 
total susceptibility. 



4-1-2. Scaling properties and Kondo temperature As in the FEKM, one might expect 
the impurity susceptibility to be a scaling function of T /Tk up to logarithmic corrections 
coming from the bulk marginal operator g. More precisely, weak and strong coupling 
regimes are characterized by the following behaviors: 

In the weak coupling regime, the lowest order perturbative expansion 

Ximp = ^(l-l/ln(T/T,,) + ...), (4.3) 

has been used to start building the scaling plot for Tk x Ximp, as shown in Fig. 8 where 
the data of Fig. 7 have been used to compute Ximp- Step by step, the temperature axis 
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has been rescaled with the Kondo temperature hke T T /Tk in order to produce 
the best collapse of the data onto a curve which reproduces the expected behavior for 
the usual Kondo problem (see Ref. [2] for instance). Ignoring scaling violations arising 
from the bulk marginal coupling (which we expect to be small), one expects the scaling 
plot shown in Fig. 8 to give a rather good estimate of the Kondo temperature Tk for 
the Heisenberg spin chain. The resulting estimates for Tk (and thus = Vs/Tk with 
Vg = 7^/2 for the Heisenberg chain) are listed in the table 3, and plotted in the inset of 
Fig. 8. One sees the Kondo energy scale Tk going to zero exponentially when — > 0. 
We have successfully compared these estimates to the FZ form [8] 

Tk oc exp(-7ry/l/Jj,-l), (4.4) 

which in the limit of small Kondo coupling Jj^ is equivalent to the exponential square- 
root form Eq. (2.37) derived in the previous section from the RG. 




I I I I I I I I I I I 

0.01 0.1 1 10 100 1000 

Figure 8. Impurity susceptibility Ximp obtained by Quantum Monte Carlo 
calculations of the SCKM (J2 = 0) with chains of length L ~ 512 coupled to an 
impurity spin with various Kondo couplings J'j^ indicated on the plot. The collapse 
for Ximp X Tk was obtained starting from the weak coupling regime with the curve 
Ximp(T' > Tk) = 0.25/T(l - l/\n{T/TK) (red line on the right) serving as a support 
for the rest of the collapse. The strong couphng limit Ximp{T <C Tk) = 0.25/Tk is 
shown by the black horizontal dashed line on the right. Inset: Kondo temperature 
estimates obtained from the collapse on the main panel plotted versus the coupling 
J'j^. The solid curve is the expression Eq. (4.4) with a prefactor ~ 5. 



Kondo Effect in Spin Chains 



22 



J'k 
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Table 3. Kondo temperature Tk and the associated Kondo length scale ^a' = tt/ (2Tr-) 
for the SCKM ( J2 = 0), Eq. (1.2), with a Kondo coupling J'j^. The values are estimated 
from the impurity susceptibility data (see Fig. 8). 



One technical remark has to be mentioned about the QMC data for Ximp- At 
the strong couphng fixed point it is quite difficuh to simulate the system at very low 
temperature T ^ Tk and in the scaling limit T ^ Vs/ L at the same time. Finally, as 
briefiy mentioned above, we expect that the bulk marginal interaction will lead to small 
corrections to the susceptibility. 



4-2. Bethe ansatz solution 

We start from an open Heisenberg chain with a boundary spin impurity Eqs. (1.2) and 
an external magnetic field B = Bel 

L-l L 

n = Y,S,- S,+i + J'Ap -Si-Y^ BSt. (4.5) 

■t=l i=l 

This model has been shown to be integrable [8] and allows a Bethe Ansatz solution 
which, in the case where the impurity spin Simp — 1/2, reads 

M 

[ei(Aj)]2^ei(Aj- + Ax)ei(Aj- - A^) = JJ ^2(Aj - Afc)e2(Aj- + A^). 

k=l,k^j 

(4.6) 

M is the number of down spins, the function e„(x) = ^I'l™^^ and the impurity coupling 
constant Jj^ is related to by 

J'k = -i^r- (4-7) 



1 + A2 



K 



Note that the \k defined here is not the same as the \k defined in Sec. II. It is more 
convenient to re- write Eq. (4.6) as 



A- -I- ^ 
A . _ i 



2L 



Aj + Xx + I 



Aj - Xk + ^ 
Aj - Xk - ^ 

M 

n 



Aj + Ak + i 
Aj + Ak-i 



Aj -Ak + i 
Aj - Ak-i 



Then we get 



2L9{2Aj) + 9{2Aj + 2Aa') + 0{2Aj - 2Xk) 



M 



J2 d{A, + A,,) + d{Aj - Ak) - irl 
k=i,kj^j 



J' 



(4.9) 
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where d{x) = arctan (x), the Ij are all integers and j = 1, M. The total energy in an 
external magnetic field B can be expressed [8] 

E(M) = ^ + :^--y ^.'^ , -(^^±1-m]b. (4.10) 
^ ^ 4 4 2 ^ A? + 1/4 V 2 / ^ ^ 



4-3. Singlet-triplet gap 

For the J2 = spin chain with open boundary conditions, the result Eq. (3.13) is 

modified, to first order in the marginally irrelevant bulk coupling constant, g{L) to [18]: 

77-7; 

Ast{L) = ^[1 - g{L)] (4.11) 

We first compute Ast{L) in the clean case (Jj^ = 1) using the Bethe Ansatz solution 
up to L = 10^ to extract the marginal coupling g{L) from Eq. (4.11). A very good 
agreement (for L ^ 100) with g{L) from Eq. (2.23) is observed using Li ^ 0.92. 

For J/^ 7^ 1 we expect Eq. (4.11) to continue to hold in the strong coupling limit of 
the SCKM, ^ L, but we may also include the Fermi liquid correction of Eq. (3.14): 



Ast{L > ^ir) = 



9{L)- 



(4.12) 



2L 

Note that this result is to first order in g{L) and first order in ^k/L only. It should be 
corrected by a Taylor series in these 2 coupling constants. We also emphasize that the 
g{L) ~ l/ln(L) term is the dominant correction in this expression. Hence, in order to 
verify the presence of the term proportional to it is natural to study the following 




Figure 9. Bethc Ansatz results for the SCKM (J2 = 0). Data are shown for 
the difference &{Jx = 1) ^ ©(^k) [Eq- (4.f3)] computed exactly for system sizes 
20 < L < 4000 with an impurity coupling J'j^ = 0.5. The solid line represents a fit to 
the data of the form Eq. (4.13) with ^a' = 6.75, a = -10.21, 5 = -109.39. 
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quantity: 

:^A'sAJ'k) = [Ast{J'k = 1) - Ast{J'k)] 

^e(j^ = i)-e(j;,) 

2Z + LM^ + L^' ^^-^^^ 

where we have included the leading terms in the expected Taylor series in the 2 coupling 
constants. In particular, one might expect a term of the form ~ {irvg/ L)g{L)C,K/ L ~ 
l/(L^ln(L)) to be present in the finite-size corrections to Ast{L). The analytical 
calculation of the coefficient of such a term would be challenging since it would involve 
diagrams with 6 current operators. In the opposite limit L <^ C,k can calculate 
Ast{L) in perturbation theory in A^-, Eq. (3.15). At weak coupling (L ^ C^k) there 
is no correction of 0{g). This follows because both singlet and triplet states arise from 
the same state (with 5* = 1/2) of the chain of L spins. 

Exact results available form the Bethe ansatz solution of finite systems allow us 
to verify the presence of the ti^k/{'^L) term in the singlet-triplet gap, Eqs. (4.12), 
(4.13). Our results for 6(Jj^ = 1) — 6(Jj^) are shown in Fig. 9 for J'^ = 0.5, clearly 
demonstrating the presence of a linear term ~ {^k ~^k{J'k = 1))/-^- If the higher order 
corrections described by Eq. (4.13) are included a very precise absolute estimate for 
— ^k{Jk = 1) can be obtained. We expect ^k{J'k = 1) to be 0{1) and therefore to 
be negligible compared to C,k for large enough C,k- Fits to the form given by Eq. (4.13) 




Figure 10. Bcthe Ansatz results for the SCKM (J2 = 0). Data arc shown for the 
singlet-triplet gap Ast versus L/{£^k — £.k{Jk = 1)) computed exactly for various 
system sizes 20 < L < 4000 and impurity couplings J'j^ = 0.3, . . . , 0.9 as indicated 
in the plot by the different symbols, ~ £.k{Jk — 1) is obtained from fitting to 
Eq. (4.13) as listed in table 4. Due the the presence of the marginal coupling, the 
data clearly do not scale with i/^A' as was the case at the critical point J| as shown 
in Fig. 5. 
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J' 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



U-^k{J'k = ^) 38.12 14.32 6.75 3.53 1.90 0.96 0.38 



Table 4. Bcthc Ansatz estimates for the Kondo length scale — £,k{J'k ~ 1) of the 
SCKM Hamiltonian Eq. (1.2), J2 = 0, with a Kondo coupling J^. 



are extremely good down to Jj^ = 0.3 and the obtained values of C,k — ^k{J'k — 1) 
listed in table 4. For smaller values of the diverging prohibits reliable fits for the 
system sizes where we can numerically solve the Bethe ansatz equations. 

As we showed in Fig. 5, the quantity 6 = L/ {ttVs)Ast{L, J'j^) scales with L/C,k 
in the FEKM and also for the SCKM at J| where the bulk marginal coupling, g, is 
zero. At J2 = 0, is non-zero and is in fact the dominant correction term in singlet- 
triplet gap, Eq. (4.12). Hence, for J2 = we no longer expect scaling of with L/^k, 
but instead, should approach 1 — g{L) for large L. This is demonstrated in Fig. 10 
where exact Bethe ansatz results for for the SCKM at J2 = are plotted versus 
L/ {^K — ^k{Jk = 1) with C,K — ^k{Jk = 1) obtained from fitting to Eq. (4.13) as hsted 
in table 4. Clearly the data do not scale, even for sizable values of — ^k{J'k — 1) 
where ^k{J'k = 1) can be neglected, but approach the value 1 — g{L) plotted as the 
dashed line with g{L) from Eq. (2.23) using the previously determined Li ~ 0.92. Data 
are shown for 20 < L < 4000. From Eq. (4.12) it would be tempting to assume that 
the quantity Q{Jk = 1) — ^{Jr) should scale with L/^k for the SCKM at J2 = 0, 
however, higher order terms in g{L) not included in Eq. (4.12) do not cancel, corrupting 
the scaling. Likewise, we also do not expect to scale with L/C,k in the weak coupling 
regime L <^ since, even though there is no correction of 0{g), since \k{L), Eq. (2.9), 
is not a pure scaling function of ^k/ unlike in the FEKM. 

4.4- Finite size scaling of the ground state energy in the SCKM with J2 = 

In section 4.4, Eq. (3.19), we derived the finite size corrections to the ground-state 
energy neglecting the bulk marginal operator, g. We now discuss the more general case 
where g > and logarithmic corrections arising from this operator are present. For the 
uniform ( = 1) chain with open boundary conditions these corrections are known [18]. 
In addition, as outlined in Appendix A, it is also possible to derive the universal term 
in the ground-state energy arising from the coupling between g and C,k'- 



Combining the various results, we find for general J2 < J| and Jj^ 7^ 1 with L ^k'- 



6Egs 



32L2 



(4.14) 



EGsiL,J'K) = eoiJ'K) + eiL 



■KVs 



1 



9g 



2AL 



2 




(4.15) 



Here, as L ^ 00, (7 ~ 1/lnL. 
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We have attempted to verify the correction to the ground-state energy, Eq. (4.14), 
for non-zero g through the use of the Bethe ansatz solution for the SCKM at J2 = [8]. 
Ahhough a logarithmic term of the correct order can clearly be identified, the sign 
appears to be the opposite of Eq. (4.14). This may be due to significant higher order 
corrections that could change the effective sign observed in the numerical results. 



4-5. Bethe ansatz results for the magnetization, M{H) 

FZ gave results on the magnetization for the SCKM at zero temperature at infinite 
length using the Bethe Ansatz. They found [8, 22] that, at large L, the magnetization 
had the form: 

M{L) = m,L + (M|,^, + M^^) + 0(1/1). (4.16) 

Here mo is the bulk magnetization per unit length, which is linear in H (up to logarithmic 
corrections) at H <^ J. For H <^ J, FZ found: 



M. 



FZ wrote M 



"'^^^ 4 [\n{Ho/H) 
in terms of a field 



l ln[(l/2)ln(/7o/g)] 
2 \n\Ho/H) 



imp 

Here: 

and \k is defined by: 



1 + A2 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



K 



They identified with a Kondo temperature, Tk- At Jj^ <^ 1 this has the exponential 
square- root form derived in section 2.4 from the RG: 



Tk = i/oe""vV^. 



(4.21) 



At H <^ Tk FZ found that M, 
Tk <^ H <^ J, they found: 

1 1 



imp vanishes (for an S=l/2 impurity). However, for 
, \n[{l/2)ln{H/TK)r 



1 



HH/Tk] 

1 



\ti\H/Tk) 

1 ln[(l/2) HHI/HTk)] 

2 \n\Hl/HTK) 



(4.22) 



\Yi{Hl/HTK) 
{Tk<^H <€. J). 

For the model with J^^ = 1, Furusaki and Hikihara [21] observed that M^^^^ (which 
was obtained earlier by Bethe ansatz for the = 1 case in related models [23, 24]) can 
be derived from renormalization group improved first order perturbation theory in g. It 
is simply {l/A)g{H)^ up to a higher order term (at small H) of 0[1/ \y?{Hq/ H)]. We 
should associate terms {l/S)g{H) with each boundary. 
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It is interesting to consider the "impurity magnetization" M.imp as it is usually 
defined in the theory of the Kondo effect. This is the additional magnetization resulting 
from adding the impurity. For the spin chain this is naturally defined as: 

M^mp = lim [M{Jk, L + 1)- M{Jk = 1, L)]. (4.23) 

Ignoring a term mQ{H), the bulk magnetization per unit length, which is smaller by a 
factor of (H/J) than the other contributions, we see that: 

M^mp = M^p. (4.24) 

The first line in Aiimp in Eq- (4.22) is precisely the standard result for the FEKM, at 
Tk <C if -C -D, up to corrections of 0[1/ \'d.^{H/Tk)]- (See, for example, [2].) It is, of 
course, a scaling function of H/Tk- The second line in M.imp in Eq. (4.22) is rather 
complicated. Using: 

HHTk/H^o) = 2 HH/Ho) + HTk/H) (4.25) 
we can write it in terms of g{H) and Xk{H). 

mmp ^ 0/ (m~l{\ (nv ^^-26) 
^I9[H) + 1/\k[H) 

This could be Taylor expanded in powers of g{H) and indeed has the form of Eq. (2.38). 
However, we find the peculiar non-analytic form of Eq. (4.26) puzzling. 

1 9^f-9 , ^427) 



2/g+l/\K 2^V2Ax 
Clearly, this expression does not have a joint Taylor expansion in g and Ax- 



5. Conclusions 

Our main results for S,k = Vg/Tx are summarized in Fig. 11 for the SCKM at J2 = 0, 
panel (a), and J2 = J|, panel (b). We have shown that the spin chain Kondo model with 
the second neighbor coupling, J2 adjusted to the critical point has the same low energy, 
long distance behavior as the usual free electron Kondo model. This leads to the usual 
exponential divergence of with 1/J^ as shown in panel (b) where the concordance 
of as obtained from scaling of A^t, Egs and RG is demonstrated. Furthermore, 
at J2 = J| many quantities display scaling with L/^k in an identical manner for the 
SCKM and FEKM (Fig. 5). For J2 < J2 ^ new kind of Kondo effect arises in the 
SCKM due to the presence of the non-zero marginal coupling g. As shown in panel (a) 
for J2 = this leads to a more slowly diverging ~ exp(-^l/ Jj^ — 1) in agreement 
with the analytical Bethe ansatz result [8]. The non-zero marginal coupling g > Q 
destroys scaling with ^k/L and reliable estimates for can only be obtained from a 
numerical determination of finite-size corrections. However, the strong coupling fixed 
point remains that of a completely screened impurity and the usual Kondo physics is 
still present in this case albeit obscured by significant logarithmic corrections arising 
from the bulk marginal operator, g. Hence, we expect physical quantities to behave 




Figure 11. (a) Results for the nearest neighbor SCKM, J2 = 0. £,k = Vs/Tk as 
obtained from the impurity susceptibihty, Ximp (table 3) shown with obtained 
from scaling of Ast (table 4). Both quantities are seen to agree with v/Tk from FZ, 
Eq. (4.21). (b) Results for the next-nearest neighbor SCKM with J2 = J|. (,k are 
shown as obtained from scaling of Ast and from the ground-state energy. Agreement 
is observed with the RG result, Eq. (2.10). 



the same way as for the FEKM with corrections which only vanish logarithmically with 
decreasing energy or increasing system size: 

/(L, J'k) = h{L/U) + g{L)h{L/iK) + g\L)h{L/iK) + • • • , (5.1) 

where fQ{L/C,K) is the universal scaling function occurring for the FEKM. 

Several open problems pertaining to the Kondo effect as it occurs in the SCKM at 
= ((7 > 0) remain. First, as briefly mentioned in section 4.1 we expect that the 
bulk marginal interaction will give rise to small corrections to the susceptibility, precise 
analytical or numerical results for such corrections would be desirable. Secondly, in 
section 4.4 we derived an analytical result for the coefficient of the {nVs/ L)g{L)^K/ L 
term in the ground-state energy, Eq. (4.14), lending support to the proposed form of 
the scaling corrections Eq. (2.38). However, even though our numerical solution of the 
Bethe ansatz equations show a logarithmic term of the correct order the sign appears 
to be the opposite of that of Eq. (4.14). A possible explanation for this would be that 
higher order are significant and change the effective sign seen in the numerics. Further 
analytical insight to the higher order terms or drastically improved numerical results 
would be valuable to resolve this questions. Thirdly, additional support for the proposed 
scaling corrections, Eq. (2.38), in the form of an analytical calculation of the coefficient 
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of the {tiv / L)g{L)^K / L term in Ast would also be of considerable interest. We have 
so far not been able to derive this coefficient in a simple manner. As pointed out in 
section 4.5, we find the peculiar non-analytic form of Eq. (4.26) puzzling since it does 
not have a joint Taylor expansion in g and Xk-, an issue we hope will be resolved in the 
near future. Finally, we hope that the spin chain Kondo models studied in this paper 
will help experimentalists for a better understanding of impurity effects in spin chains 
materials. 
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Appendix A. Derivation of g{L)$^K/L^ term in ground state energy 

We follow closely [18], which we refer to as AQ, but note the change in notation: 

gAQ = {V3/A7T)g, Jaq = {II2ti)J. (A.l) 

First let us correct (2.21) of AQ. This was meant to be the connected part but the 
subscript "connected", as in (2.18) was missing. More substantially, the connected 
Green's function should have two other terms, missing from (2.21). The correct result, 
with J normalized as we do is: 



c 



ZiZ2{t - Zi){t - Z2) 



^ab ^ab 

(We dropped the L subscripts but all currents are left-moving.) These last 2 terms, 
missing from (2.21) of AQ, do not contribute to the g^ / L term in Egs although they 
do lead to an additional non-universal g^ term without the 1/L factor. Therefore, they 
make no change in the conclusions of AQ. 

Now consider including the Fermi liquid interaction as well as the bulk marginal 
interaction: 



L 



H = Ho- ^riO) - f- / dxJix) ■ Ji-x). (A.3) 



For infinite L, standard perturbation theory to first order in g and C,k gives: 

6 2n 



6Egs = -^^ I dr I dx{f{0)J{x, t) ■ J{-x, r))c. (A.4) 
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From Eq. (A. 2) setting r — ^ 0, 21 — * r + 2X, 2:2 ^ t — ix, and summing over a = b, this 
is: 

SEgs = -^-f^ r dr r dx I --. (A.5) 

6 27r Jq 2{t + txy{T - txy 

Note that all 3 terms in Eq. (A. 2) contribute with the 2 last terms giving 1/2 the 
contribution of the first term. Also note that the integrand is an even function of x so 
we can extend the region of integration to — cxo < a; < 00. Doing the x- integral gives: 

6 27r 8 |r|3 ^ ^ 

This integral is divergent. We cut it off, |r| > tq, giving: 

SEo. ^ - (A.7, 

This is a non-universal term with no L-dependence. To get the L-dependent term we 
need the generalization of Eq. (A. 2) to finite L. This is obtained by: 

T±ix ^ {2L/n) sinh[(r ± ix)Ti/2L\. (A.8) 

the X and r integrals (with the cut off on the r integral) can again be done exactly. 
Again the integrand is an even function of x so we can extend it to —L < x < L and 
then change variables to z = e*'^^'/^. The resulting contour integral is elementary, with 
a double pole inside the contour giving: 

32L'^ J-oo I smh(7rr/L)|-^ 
Doing the cut-off r integral exactly gives: 

^^^^ = "32[(L/vr)^nh(7rro/L)]2- (^'^^^ 
Now we Taylor expand to second order in tq/L giving: 

oEgs 



32 

So the desired universal term is: 



LTo^ 3L2 + • • ; 



(A.ll) 



Reinserting the factor of Vs which was set to 1 gives the 0{g/L'^) term in Eq. (4.15). 



Appendix B. Extracting the l/L^ Term in Eqs 

In this appendix we describe a form of Richardson extrapolation that is useful for 
isolating the 1/L^ term in the ground-state energy needed to extract in Eq- (3.19). 
We assume that the ground-state energy is of the form, Eq. (3.19): 

eW(^) = i(L) + 5c°(L) + Cc^iL) + DcliL) + OiL~''), (B.l) 
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with c\{L) = L,cO(L) = l,cO(L) = 1/L and cO(L) = l/L^. We're interested in 
determining the coefficient D. Furthermore, we assume that E{L) is known for a 
sequence of L. It is then easiest to proceed using a form Richardson extrapolation 
by systematically removing the A,B,C terms. We eliminate the term proportional to 
A by constructing a new series: 

= - + 2)^0^^^ (B.2) 

Naturally, E^^^ must be of the form: 

E^^\L) = BcliL) + CcliL) + Dc\{L) + 0{L-^), (B.3) 

with 

c!(^) = c°(L)-c°(L + 2)-;^^^^^, ^ = 0,1,2. (B.4) 

Now we can eliminate the term proportional to B by transforming: 

E^^\L) = E('\L) - E('\L + 2)^f^, (B.5) 

which we can write as: 

= CcliL) + DcliL^) + 0{L-^) 
cKL) = c\{L) - c]{L + 2)^^ ^ = 1^2. 



Repeating once more we obtain E^^'^ (L) of the form: 
E'^^\L) = DcliL^) + 0{L~^) 
cliL) =cl{L)-cl{L + 2)^^^. 

From which the constant D can be extracted by plotting E^^\L) / C2{L) . 



(B.6) 



(B.7) 
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